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. Abstract 

: 

The mass formulae for finite lumps of strange quark matter with u, d and s 

o ■ 

H ■ 

' quarks, and non-strange quark matter consisting of u and d quarks are derived in 

: 

J> . a non-relativistic potential model. The finite-size effects comprising the surface, 

' curvature and even, the Gauss curvature were consistently obtained, which shows 

a converging trend. It is found that there is a possibility for the formation of 
metastable strangelets of large mass. The model predicts low charge to mass ratio 
as the characteristic signature of strange matter in agreement with the relativistic 
studies. This study also yields an independent estimate for the bag energy density 
which is in agreement with the M.I.T bag model value. 



I. INTRODUCTION 



The possibility of strange quark matter to be the ground state of QCD and the probable 
stability of such matter in finite lumps called as strangelets, speculated by many authors[l- 
5] has generated a spurt of activity during the last few years. The general framework 
adopted for such studies has been the relativistic Fermi gas model with relevant QCD 
parameters like the bag constant B, strong coupling constant Oc and the strange quark 
mass rris. Due to large uncertainty in these parameters, it has not yet been possible to 
arrive at a conclusive result. 

The calculations of Chin and Kerman and Freedman and Mclerran [^], which 
includes the lowest order quark-quark interaction shows that infinite quark matter even 
with strangeness is unstable against the nucleon decay. Regarding the probable stability 
of a finite lump of quark matter, Farhi and Jaffe were the first to consider the mass 
of the strangelets by adding the surface energy and the Coulomb energy terms to the 
volume energy, by initially treating the surface tension coefficient as a parameter. Later 
Berger and Jaffe |^ developed a formula to calculate the surface tension coefficient. They 
showed that surface energy is zero for massless case, and to be about 100 MeV for the 
massive case [m„ = m^^ = 0, ~ ISOMeV"]. In this work, Berger and Jaffe have dropped 
the quark-quark interaction with the supposition that it can be absorbed into the bag 
constant parameter. Thus, they have developed a mass formula assuming a system of 
non-interacting quarks confined in a finite-size cavity for which they included the surface 
and the Coulomb effects. 



Recently, following Mardor and Svetitsky [0 who have pointed out that the curvature 
energy might play a decisive role for the stability of strangelets, Madsen has obtained 
an estimate for the curvature energy coefficient to be about 435 MeV for the massless 
case. This has dramatic consequences for the stability of low mass strangelets. In view of 
its importance, the curvature coefficient should also be determined for the massive case, 
which has not yet been derived. The indication of predominance of the curvature over 
the surface term is already intriguing. Therefore, it is highly desirable to develop a mass 
formula with surface, curvature and Gauss curvature terms correctly calculated showing 
convergence trend. During the last two decades, extensive application of non-relativistic 
potential models to the description of ground-state, excited-state and spectroscopic 
properties of baryons and mesons have been amazingly successful. In view of this, and 
due to the momentous nature of the problem of stability of strangelets and its far reaching 
consequences, it is worthwhile to attempt at a development of a mass formula based on 
a non-relativistic approach. In this paper, we report on such an attempt, which has the 
advantage of reliable estimation of the various coefficients with practically no parameter 
except the baryon number density. 

In Sec. II, we present the details of our non-relativistic formalism used to derive the 
mass formulae for multiquark droplets including the finite-size effects. We then discuss the 
results obtained in our model in Sec. III. Finally, we summarize and present our conclusions 
in Sec. IV. 

II. FORMALISM 



In this work, similar to Berger and Jaffe relativistic study, we consider a system of 
finite number of non-interacting particles confined in a box. This can be represented by 
an infinite square well potential for which the correct enumeration of quantum states has 



been done by Hill and Wheeler who obtained the density of states for such a system 



as 

dN = ^k^dk - —kdk + —dk (1) 

where g is the degeneracy factor. For a spherical box of radius R; V = |7ri?^, 5* = 47ri?^ 
and L = 27cR. It is then straight forward to obtain an expression for the total number of 
particles as 

iv=r^* = 44-^4 + ^*F (2) 
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and for the total kinetic energy Ekm as 

fkFh^edN,, 
Ekin = I —dk = — 
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where kp is the Fermi momentum and m is the mass of the particles. Considering a 
system of one kind of quarks, we have g = Q. We then express the Fermi momentum kp 
in terms of the radius R, by inverting Eq.(2) and retaining terms upto 0{N~^) as 

kp ^ {7r'py/'[l + Ci(7rV)-'/' + C^irc'p)-'/' + dn'p)-']'/' (4) 

where Ci = 3nS/{8V), C2 = 2Cf/3 - 3nL/{4V), C3 = Cf/3 - Ci3'kL/{AV), and the 
total particle number density p = N/V. Then, substituting the above expression for kp 
in Eq.(3), and grouping all the terms with same powers of N, we arrive at the following 
expression for Ekin upto the order of 



The coefficient of N, N'^^^, N^^'^ and are termed as volume, surface, curvature and 
Gauss curvature contributions respectively to the total kinetic energy. 

Now, in order to obtain an expression for the total mass of a multiquark system, the 
following important effects are to be included. 

i) For the system to be stable at finite density, pressure has to be balanced at the 
surface. To ensure this, we have to add an extra term BV which is quite analogous to 



the MIT bag term |11 . 



ii) The center of mass correction for a many-particle system can be taken into account 
in an approximate way by multiplying the factor /2m in the kinetic energy expression 



by a factor {^^jf-)-, as customarily done in the non-relativistic Hartee-Fock calculations 
[0- 

iii) In addition to the finite-size effects, the Coulomb effect has to be included . For a 
spherical system of radius R and charge Ze, the Coulomb energy is given by Ecoui = ^f^f"- 

iv) Lastly, the rest mass energy of all the particles (= mN) confined in the spherical 
box must be included. 

Thus, including all the above effects, the total mass of a particle system is given in 
a compact form as 

M{N) = Ekin + BV + Ecoui + rnN. (6) 

We then eliminate the bag energy density B by using the stability condition, which is 
nothing but the pressure balance equation 

d{M/N) 



dp 
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0. (7) 



where poq is the equihbrium quark number density. This gives, 
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Then, in the bulk hmit, neglecting the Coulomb effects, we can rewrite B as 



Thus, we arrive at a formula for the bag energy density B which can be computed using 
the mass of the constituent quarks and the number density of the system. 

It is interesting to recall here that the usual bag models of hadrons employs a bag 
constant B, which is considered an universal constant of nature. It is being widely used in 
the study of hadron properties and the phase transition of hadronic matter to quark-gluon- 
plasma. Normally, the value of B has been determined from the ground state properties 
of hadrons, which hes between 58 MeVfm~^ [B^l^ = U5MeV] and 86MeVfm-^ [S^^ = 
IGOMeV]. 

On the nonrelativistic front, potential model studies have also been very successful in 
describing the baryon properties. Therefore, the above formula for B given by Eq.(9) is 
most likely to give an independent realistic estimate of the same. We then express the 
quark number density Pog in terms of the baryon number density Po as Po = PoqI'i and fix 
the mass of the w-quark to be 300 MeV which is normally used in non-relativistic models 



[|T^.We take the value of po to be about 1.3pn„c, where the nuclear matter density pnuc 
is given by 0.153 fm~^ WM. Using these values, we obtain our value of B to be about 



83 MeV fm ^. Amazingly, this value lies within the range of relativistic values of about 
58 MeVfm-^ and 86 MeVfm-^. This result gives us confidence in the reliability of our 



model for the study of stability of multiquark droplets. Further, the above formula given 
by Eq.(9), clearly brings out the dependence of the bag energy density B on the baryon 
number density and the quark mass. 

Finally, we now obtain the following mass formula for a system of N particles, each of 
mass m and charge q, by substituting the value of B given by Eq.(8) in Eq.(6). 



The mass formula given by Eq.(lO) is for a system consisting of one kind of quarks. It 
can be generalised to obtain the mass of a system consisting of quarks of more than one 
flavour. Here, we apply it to a multiquark droplet consisting of u and d quarks, hereafter 
referred to as "udilet" . For a given number of u and d quarks and respectively, we 
define an asymmetry parameter S = {Nd — Nu)/Nq, where Ng = Nu + Nd- So, the quark 
density for each flavour can be expressed in terms of S and poq as 



with = Ng/V . The total kinetic energy for a two flavour system is obtained by 
summing over the respective contributions, given by Eq.(5) for u and d quarks. Then, one 
can obtain the total mass of an udilet in terms of the baryon number density Po[= Poq/'i] 
using Eqs.(lO) and (11) as 



M{N) = -Ekin + -Ecoui + mN 



(10) 



A. Two flavour system. 
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where we have used rrid = rriu and A = By substituting for p^, and the radius 

R — ToA^I^ with ^nr^po — 1, we thus arrive at the Bethe-Weizsacker hke mass formula 
for the udilets upto the order of and S^. 

= a,2 + as2A-'/' + a,2A-'/' + 0^2^-^ + a^Z^A''/' + 3m„ (13) 

where a^2,as2 etc are defined as a^2 = civi^ + |^^)) 0fs2 = ^^(l + §5^), ac2 = Oc, ^02 = 
ao(l — |(^^), and the various coefficients are given by, 

5 971^/3 



3 2m^ 

5 /3,c;/Q,3 



2/3 



-)5/3(-)4/32^2^2/3 



3 2mu r '2 

3 2;;^[^[^-2-]-^]^^' 

4.3e2, 



B. Three flavour system 
Similar to the two flavour system discussed above , we now consider a muhiquark 
droplet consisting of u, d and s quarks usually referred to as strangelet. For this case 
of three flavour system, we need to deflne two asymmetry parameters namely, S^d ~ 
(Nd - Ny,) /Ng and 5us = (Ng - Ny,) /Ng, where Nq^Nu + Nd + Ng. We can then express 
the quark density of each flavour in terms of dud, Sus and poq as 

Pu = - Sus - Sud) 

Pd = - 5us + 25ud) 

= P^(i + 26us-Sud). (15) 

The total kinetic energy for this system is obtained by summing over the respective 
contributions given by Eq.(5) for m, d and s quarks. Using Eqs.(lO) and (15), we thus 
arrive at the following mass formula for strangelets in terms of the baryon number A 
and its density Po retaining, terms upto 0{A~^) and second order in the asymmetry 
parameters. 

= a,3 + assA-'/' + a^sA-^/' + a^sA-' + a^Z^A'^/' + (/„ + /,)m„ + /,m, (16) 

where the quark fractions can be expressed using Eq.(15), as follows . 

N 

f — — 1 _ A _ A 

Ju — — Oud 

fd — — ^ — ^us + '^^ud 

N 

fs = = 1 + 25,3 (17) 



The various coefficients a„3,as3 etc in Eq.(16) are given by, 



where 
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Ai = {6ud + ^usf + i'^Sud - Susf 

Al = {26us-Kdf. (19) 

III. RESULTS AND DISCUSSIONS. 

The mass formulae given by Eqs.(13) and (16) are most general, which can be used 
for any arbitrary combination of u,d and s quarks. To calculate the various coefficients 
like volume, surface etc , we have chosen the mass of u and d quarks to be 300 MeV 
which lies within the range of values normally used in the non-relativistic potential model 
calculations |]T3 of baryon spectroscopy. We take the mass of s quark to be 400 MeV, 



which is quite reasonable. The baryon number density po is assumed to be about 1.3 



times the nuclear matter density pnuc = 0.153/m~'^. For the sake of convenience, we 
calculate these coefficients by imposing charge neutrality condition, which for udilets 
imphes A^^ = 2A^„, 5 = | and for strangelets Nu = Nd = Ns, Sud = 0,6us = 0. The 
value of the coefficients so computed for udilets and strangelets are given in Table I. It is 
satisfying to note that the contribution of the various terms successively decreases leading 
to a faster convergence of M{A)/A. It is interesting to note that the Gauss curvature 
term has negative sign and so has the opposite effect on the mass compared to surface 
and curvature terms. This is in fact is a welcome feature. It must be realised that these 
kind of mass formulae based on Fermi gas approximation are not valid for very small 
baryon numbers. Because of its general power series structure in A~^^^, the higher order 
contributions becomes progressively important for low A values. However, we feel these 



formulae are more appropriate for baryon number A > 50 as also noted by Jaffe [|15 
In Table II, we have presented the corresponding masses of udilets and strangelets for 
certain typical values of A. From symmetry considerations, strangelets are expected to 
be more stable than the corresponding udilets for a given baryon number. This is borne 
out quite well in our calculations as can be seen from Table II. The energy per baryon for 
strange quark matter is about 1275 MeV, and for non-strange quark matter consisting 
of u and d quarks , it is about 1320 MeV. Thus, both these quark matters are unstable 
against nucleon decay. So, we find in this model the nuclear matter to be the ground 
state of hadronic matter. However, this does not apriori preclude the possibility of meta- 
stable strangelet formation. Before, we embark upon the investigation of the possibility of 
strangelet formation in this model, we would like to consider the influence of the baryon 



number density po and strangeness on the relative stability of strangelets and udilets. 

In Fig.l, we illustrate the effect of Po on the relative stability between the strangelets 
and udilets, using the general mass formula for strangelets and udilets given by Eqs.(16) 
and (13). We have found the most flavour-stable isobar, by minimising the mass with 
respect to 5 for udilets, and 5ud and 5 us for strangelets. In Fig. la, the difference between 
the masses of the strangelets and udilets so obtained are plotted as a function of po/ Pnuc, 
for a typical baryon number A — 400. It can be seen that mass of the strangelets is lower 
than the corresponding udilets for all densities. The minimum difference between their 
masses is about 4 MeV. This feature is true for any baryon number A is shown in Fig. lb, 
where the surface and Coulomb effects are neglected 

In Fig. 2, we illustrate the effect of strangeness on the relative stability between the 
strangelets and utilets for the same A — 400. Using the general formula for strangelets 
given by Eq.(16), we have determined for a given Sus, the flavour-stable value of mass 
Ms^A) by minimizing the toal mass with respect to 6ud- For two chosen values of po, 
namely 1.27 Pnuc a-^d 2.14 Pnuc, the results so obtained are presented as solid curves in 
Fig. 2 as a function of fg, where fg in terms of dud and dus is given by Eq.(17). It is 
interesting to flnd that the mass shows a well deflned minimum for fg at about 0.7. Then 
to show the relative stability with respect to udilets having the same baryon number A, we 
have minimised the mass M2{A) given by Eq.(13) with respect to S. We have represented 
the most stable mass of udilets so obtained for the same two values of po as dashed lines. 
It can be seen from Fig. 2 that at the optimum value of fg, the difference between the 
mass of strangelets and udilets is maximum, which gradually decreases and tends to zero 



as fs is varied about the optimum value. Therefore, there is a maximum value of fs, 
referred to as /J", above which udilets becomes more stable than the strangelets. Thus, 
there is a 'window' in fs space, only within which strangelets are more stable than the 
corresponding udilets. This window region depends upon density and increases as Po is 
increased, which implies that hyper strange multiquark droplets are more probable at 
higher densities. This is shown in Fig. 3, where /™ is plotted as a function of Po/Pnuc for 
the same value of A, namely A — 400. It can be seen that as the density increases, /J" 
increases. The shaded portion in Fig.3 represents the region where strangelets are more 
stable relative to udilets, and the unshaded region represents the vice-versa. 

Now, we would hke to investigate the probable stabihty of strangelets in this model. 
For a given Sus, we determine the most stable mass per baryon Mv{A)/A of the strangelet 
in the bulk limit, by minimizing the mass given by Eq.(16) with respect to 5ud- In the 
absence of Coulomb interaction — N^,; fs{— Ng/A) is determined by using these values 
of 5ud and 5us in Eq.(17). Following Chin and Kerman, we plot in Fig. 4 My/ A — 
versus fs, where the masses of the nucleon, lambda, cascade and omega are shown as solid 
dots. It can be seen that strangelets are unstable against nucleon and lambda decay even 
at the optimum value of fs, which is about 0.7. As the optimum value of fs is less than 
unity, the number of s quarks is less than the number of u and d quarks. Therefore, the 
excess of the u and d quarks over that of s quarks, will occupy the higher momentum 
states than the s quarks. These u and d quarks will then form nucleons and strongly 
decay, analogous to the nucleon dripping from a highly radioactive nucleus. At this stage, 
as the mass of the lambda is higher than that of the nucleon, we expect the nucleon decay 



to be more probable than the lambda particle. Thus, the system is likely to be driven 
towards a state of greater strangeness, where it is stable against cascade and omega decay, 
as per our calculations. Under these conditions, strangelets can be formed in a metastable 
state. It can only decay by weak leptonic decay process. This possibility gives us some 
hope of detecting strangelets, but only of very large mass. 

In the present model, we have the density as the only free parameter, because the 
masses of u, d and s quarks are constrained by the non-relativistic potential model studies 
of baryon spectroscopy. We find that with the increase of density, the masses of the quarks 
have to be lowered in order to make the system stable against omega emission. For a given 
density, the masses of these quarks can be varied within the acceptable limits to ensure 
the same. In Fig. 5, fixing the density at about 1.3 pnuc, we have shown the range of m„ 
and TUs — niu for which the omega baryon is bound in our model. It is interesting to note 
that the mass of s quark comes about to be order of 100 MeV higher than the u and d 
quark masses, which is quite reasonable. Thus, we find that there is a sizable range of 
these parameters in which metastable strangelet formation is allowed. 

As found by Farhi and Jaffe in the relativistic study of strange matter, a possible 
signature for strange matter is its characteristic low charge to mass ratio, unlike in the 
case of ordinary nuclei, where Z ~ A/2. We would like to investigate whether our model 
would corroborate their finding. So, for the case of a strangelet of baryon number A, we 
obtain the flavour-stable value of the charge number given by 

z; = -A[6:, + 6*J 



where 5*^ and 6^^ are determined by simultaneously the two equations, 




us 
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dMs/A 
dSud 







using Eq.(16). 

Similarly, the flavour-stable value of the charge number Z2 for the corresponding udilet 
is obtained using 



where, the mass of the utilets M2{A) is given by Eq.(13). 

In Table III, we have presented the values of Z2 and obtained in the above calcula- 
tion for various baryon number A. For the sake of comparison we have also given the value 
Z — A/2, which is appropriate for nuclei of large A. It can be seen that the ^min for the 
udilets {Z2) is somewhat lower than the corresponding nuclei, which becomes substantial 
with the increase of A. However, for the strangelets the Z^^^ value (Z3) is strikingly low. 
It is indeed satisfying to find that the present non-relativistic study arrives at a similar 
conclusion as that of the relativistic study of Farhi and Jaffe. Thus, the low charge to 
mass ratio as a strong signature of strangelets is being reinforced by our present study. 

In summarising our results and discussions, we find that the present non-relativistic 
model has surprisingly succeeded in describing most of the features of strangelets as 



^2 = ^[1-35*]. 



with 5* determined from the condition 



dM2/A 
d5 



5* 



= 



found in the relativistic studies. This success and the rehabihty of our model can be 
solely attributed to the sound foundation based on the following two elements, due to a 
which straight forward calculation was possible: 

(i) . The confinement of quarks is best represented by the infinite square well potential. 

(ii) . For this potential, the exact enumeration of quantum states was given by Hill 
and Wheeler, which has been used. 



IV. CONCLUSION 

In conclusion, we have developed the mass formulae for finite lump of quark matter 
for systems with u and d quarks, and also for systems with u , d and s quarks, in a 
non-rclativistic model. The complete and consistent mass formula including all finite-size 
effects comprising the surface, curvature, and even the Gauss curvature has been derived 
for both the strangelets and udilets. For the first time, we have obtained an estimate of 
the Gauss curvature term. Further, this model study provides an independent estimate 
of the bag energy density B, which we find to be in agreement with the M.I.T bag model 
value. 

We have shown that for a given baryon density there exists an upper bound on the 
strangeness below which strangelets are more stable than the udilets. This upper limit 
on the number of s quarks increases as density increases, which implies that hyperstrange 
strangelets are more probable at higher densities. 

Finally, our calculation shows that there is a finite possibility for the formation of 
metastable strangelets of very large mass. The present study also predicts low charge 
to mass ratio as the characteristic signature of strange matter in agreement with the 
relativistic studies. 

The success of this model is due to the use of infinite square well potential, which 
adequately describes the confinement of quarks and for which exact quantum density of 
states can be obtained. 
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FIGURE CAPTIONS 

FIG.l. The difference between the mass per baryon of strangelets and udilets is plotted 
against the baryon number density po in the units of the nuclear matter density pnuc'i (a) 
for A = 400 and (b) in the bulk limit neglecting the Coulomb effects. 
FIG. 2. The mass per baryon of strangelets relative to the nucleon mass is plotted against 
the s quark fraction fg as solid curves for two chosen values of baryon density po, for a 
given baryon number A — 400. The upper curve corresponds to Po — 2.14p„„c ^-nd the 
lower one to 1.27pnuc- The corresponding mass of udilets relative to the nucleon mass is 
shown as dashed lines. 

FIG. 3. The maximum value of the s quark fraction /J" is plotted against the baryon 
density po in units of the nuclear matter density Pnuc, for a given baryon number A = 400. 
The dotted region represents where the strangelets are more stable than the corresponding 
udilets. 

FIG. 4. The mass per baryon of the strangelets M^/A relative to the nucleon mass in 
the bulk hmit (neglecting the Coulomb effects) is plotted against the s quark fraction fg. 
The masses of the nucleon (Mjv), lambda particle (Ma.), cascade particle (Ms) and omega 
particle (M^) are shown as solid dots. 

FIG. 5. The range of the parameters m„ and nis, the mass of u and s quarks respectivley 
favouring stability of strange quark mass against omega emission. The contour plot of 
the mass per baryon My/A of the strangelets with /s = 3 (see Eq.l7) in the bulk limit 
(neglecting the Coulomb effects) is plotted against and difference between m„ and m^. 
The corresponding numerical values of My/A are shown . 



TABLE CAPTIONS 
Table I. The values of the volume, surface, curvature and Gauss curvature coefficients, 
Qy, Qs, ttc and Oo respectively obtained using the charge neutrality condition are given for 
both the udilets and strangelets. 

Table II. The mass per baryon for strangelets (Ms/A) and udilets {M2/A) calculated 
using the charge neutrality condition are given for various baryon number A. 
Table III. The flavour-stable value of the charge number for udilets {Z2) and strangelets 
(Z3) are given for various baryon number A. 



Table I 



System 




tts 








(MeV) 


(MeV) 


(MeV) 


(MeV) 


udilet 


417.53 


779.14 


510.59 


-286.77 


(udd) 










strangelet 


275.10 


608.89 


468.04 


-304.67 


(uds) 











Table II 



A 


M2/A 


M3//I 




{MeV) 


{MeV) 


50.00 


1560.90 


1468.77 


100.00 


1506.22 


1424.96 


200.00 


1464.25 


1391.38 


300.00 


1444.35 


1375.48 


500.00 


1423.22 


1358.63 


1000.00 


1400.26 


1340.36 


3000.00 


1373.91 


1319.46 


5000.00 


1364.78 


1312.25 


10000.00 


1354.76 


1304.34 



Table III 



A 


Z^A/2 


z* 


Z* 


300.00 


150.00 


103.79 


28.78 


500.00 


250.00 


152.78 


44.27 


1000.00 


500.00 


246.42 


74.38 


3000.00 


1500.00 


470.38 


147.40 


5000.00 


2500.00 


609.49 


193.17 


10000.00 


5000.00 


837.78 


268.86 


50000.00 


25000.00 


1588.38 


521.26 


100000.00 


50000.00 


2041.08 


674.98 


1000000.00 


500000.00 


4503.32 


1516.66 



